Berry phases in superconducting transitions 
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I generalize the concept of Berry's geometrical phase for quasicyclic Hamiltonians to the case in 
which the ground state evolves adiabatically to an excited state after one cycle, but returns to the 
ground state after an integer number of cycles. This allows to extend the charge Berry phase 7 C 
related to the macroscopic polarization, to many-body systems with fractional number of particles 
per site. Under certain conditions, j c and the spin Berry phase 7 S jump in n at the boundary of 
superconducting phases. In the extended Hubbard chain with on-site attraction U and nearest- 
neighbor interaction V at quarter filling, the transitions detected agree very well with exact results 
in two limits solved by the Bethe ansatz, and with previous numerical studies. In chains with spin 
SU(2) symmetry, ~/ 3 jumps when a spin gap opens. 
PACS. 03.65 - Quantum theory; quantum mechanics. 
PACS. 71.10 - General theories and computational techniques. 
PACS. 74.20-z - Theories and models of superconducting state. 



In the last decade, Berry phases have caused a great 
deal of interest in a variety of fields in physics. Applica- 
tions of this concept to condensed matter began with the 
study of Zak of the dynamics of a Bloch electron as its 
wave vector k changes adiabatically due to an external 
perturbation until it reaches k+G where G is a recipro- 
cal lattice vector Although in this case the Hamil- 
tonian which describes the evolution of the periodic part 
of the Bloch function is noncyclic, the initial and final 
points are related by a gauge transformation and a Berry 
phase can be defined for each band Q, which describes 
the center of gravity of the density of the Wannier func- 
tion, or Wyckoff position Later progress, showed 
that changes in the macroscopic polarization of a band 
insulator in an independent particle approximation, are 
proportional to changes in a Berry phase, which I call 7 C 
The formalism was generalized to the many-body 
case PJtJ extending previous derivations for quantized 
charge transport g, and was applied to study a ferro- 
electric transition in a strongly correlated model |J. 

Recently Gagliano and me have introduced the spin 
Berry phase j s , and have calculated 7 C and ~f s for the 
first time in a gapless metallic phase 10, 0]. The model 
was an extended Hubbard chain with correlated hopping, 
and due to the presence of inversion symmetry in it, 7 C 
and 7 S are quantized and can only have the values and 
7T. The topological quantum numbers 7 C /7T and 7 s /7r 
were used as order parameters to construct a phase dia- 
gram, separating the three phases of the model: charge- 
density wave (CDW), spin-density wave (SDW) and the 
metallic phase. The presence of anomalous flux quanti- 
zation jnj, and additional arguments p|, suggest that 
this phase is superconducting (S), possibly of triplet char- 
acter. This work opened the possibility to study phase 
diagrams in strongly-correlated systems, by looking at 
topological transitions, which are sharp even in systems 



of finite size. However, all calculations of 7 C and 7 S in 
many-body systems, were so far restricted to an integer 
number of particles per site. This restriction is too se- 
vere, particularly if one is interested in superconductors. 

In this Letter, I extend the definition of the Berry 
phase to the case in which it is necessary to perform more 
than one cycle in a parameter space, while the Hamilto- 
nian ends at a point related to the starting one by a 
gauge transformation, in order for the system to return 
to a state equivalent to the initial ground state. At first 
sight, this seems to contradict the adiabatic theorem, 
since there should be a crossing of energy levels before 
the first cycle ends. However, this crossing is not essen- 
tial if the matrix element of the Hamiltonian between the 
states involved vanishes, as it is the case if both states 
differ in some quantum number. While it seems diffi- 
cult to imagine a quantum-mechanical system with these 
properties, the generalization below of the charge Berry 
phase 7 C to systems with fractional number of particles 
per site is an example. Calculating 7 C numerically in a 
system of 12 sites, I obtain the phase diagram of the ex- 
tended Hubbard chain for 1/2 particles per site, negative 
U and any V. The results are compared with previous nu- 
merical ones, and with analytical ones obtained mapping 
the model in the limits \U\ > t, \V\ and |Z7|, V > t to 
Bethe ansatz solvable cases through appropriate canon- 
ical transformations. In one dimension (ID), for SU(2) 
symmetric systems, I show that the spin Berry phase 
jumps at the boundary between dominating singlet and 
triplet correlations at large distances, for any filling. 

For simplicity I restrict the discussion to ID and one 
band. Extension to the general case is straightforward. 
I consider a system of L sites and number of particles 
N with N/L = n/l, and n/l irreducible. The creation 
operators satisfy arbitrary boundary conditions for both 
spins Cj +La = e 1 *"^, and the Hamiltonian ff($|, 3>jj 
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conserves number of particles, z component of the total 
spin, and is invariant under translations. This invari- 
ance allows to define weighted irreducible representations 
of the translation group characterized by the total wave 
vector K = K + (7V T $ T + N^^/L, where N a is the 
number of particles with spin a and K — integer x2tt/L 
is one of the allowed wave vectors for periodic bound- 
ary conditions (PBC) 0|. Using the gauge transfor- 
mation ct = e~ l i® a / L c]~, the Hamiltonian is converted 
into one in which the fluxes $ CT arc distributed equally in 
each link -ff (<I>f , $4) (the hoppings acquire a phase M), 
with creation operators satisfying PBC c ! - +La = c*- a . 
Note that while 5(0,0) = B(±2w, ±2tt), H(0,0) = 
#(0,0) ^ H(±27r,±27r) §. To generalize 7c ( 7s ), the 
idea is to start from the ground state \gic(&i, $4.)) of 
-ff ($f, $4.) for wave vector JC and fluxes <& CT which mini- 
mize the energy, and follow the phase of \gx) as the $0- 
are shifted adiabatically by the same amount, with the 
same (opposite) sign, until \gx) reaches a state equiv- 
alent to the initial one (the same eigenstate of H ex- 
cept for a phase). I discuss first 7c . It can be defined 
as 7c = i I l^d^igK | -§^9k ($,$)>, or in the nu- 
merically gauge invariant form, discretizing the interval 
< $ < 2nl into M points $ l = 2-nli/M: 

7c = - lim Im]nUg- 2 {g K (%^)\g K ^ i+v ^ i+1 )) 

x (g K {®M-i , $ M _i) |e 1 ^ jc ^" \g K (0, 0)) (1) 

where the last ket is \gx (27rZ, 27rZ)) calculated using the 
gauge transformation from the eigenstate |<7_ff(0, 0)) of 
-ff(0, 0). In practice, 7c can be calculated to four digits 
accuracy using M ~ 6 adequately chosen points. 
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FIG. 1. Top: Scheme of the evolution of the ground state 
with flux $ for a non interacting tight-binding model, showing 
the doubly occupied one particle states k (solid circles) in the 
one-particle energy dispersion e(k). Bottom: Energy as a 
function of flux of the low-lying energy levels. 

In previous cases, I = 1 and the integration was re- 
stricted to one cycle of H (0 < $ < 2tt). However, 
while K is kept fixed, K evolves from K to K + 2wn/l 



in this cycle. Since when I 7^ 1 these K are not equiv- 
alent, the initial and final states are orthogonal and a 
Berry phase, even in its more general form p3| , can- 
not be defined for this circuit. In Fig. 1 I illustrate 
the evolution of the ground state with <!> for L = 12, 
n/l = 1/2 for the case of a non-interacting nearest- 
neighbor (NN) tight-binding Hamiltonian (Eq. (||) for 
U = V = 0). The minimum energy corresponds to 
K = $ = 0. The adiabatic continuation of this state 
in the representation of H as $ changes is always the 

Same |$o($,*)) = YI^-k/gA^I/gM > wh<3re C L is 

the Fourier transform of c^ a . In the representation of 
H , each one-particle k is shifted by $/L with respect to 
the corresponding k. It is easy to see that the energy of 
|So($,$)) is E (<f>) = £(0)cos($/L), where E(0) is the 
energy for $ = [ fT2| . As $ increases from to 2-7T, each 
k goes to the next allowed wave vector for PBC, and \go) 
evolves to an excited state with K = it. For $ near 27T, 
the ground state is = IL c -Tr/3o c -K/6* c L\°) 

with K = it. The crossing occurs for $ = it and the 
ground state energy is of course periodic with period 27r. 
This crossing does not affect the adiabatic theorem, since 
the states \go) and {g^} have different K. In contrast, |<?o) 
for <f> = ±(2ir — + ) represent orthogonal eigenstates of 
H with the same K = (see Fig. 1) and 7c is unde- 
fined (one of the factors in Eq. (1) vanishes). However, 
since both states have the same quantum numbers, they 
should hybridize when interactions are present, remov- 
ing the degeneracy at $ = ±27r. In general, for any L 
and n/l — 1/2, the above mentioned states differ in four 
particles displaced from the neighborhood of the Fermi 
point — 7r/4 to near the other Fermi point 7r/4. Umklapp 
processes of this type are generated in second and higher 
order perturbation theory in the interaction, and were 
studied previously (l4"| . For other fractional occupan- 
cies the situation is similar. The nature of the resulting 
thermodynamic phase as well as the resulting value of 7c 
depends on the detailed form of the interaction. 

I have reexamined the relation between 7c and the 
macroscopic polarization in the present case, following 
the derivation done previously for I = 1 in the many- 
body case [||-|8j . Essential for the extension is to work in 
a subspace with definite K, in such a way that for all i> 
the ground state is separated from the excited states with 
the same quantum numbers. I obtain that changes in po- 
larization AP are related to the corresponding changes 
in 7c when some parameter of the Hamiltonian varies by: 

AP = (e/27r/)A 7c (modulo e/l). (2) 

The extra factor I in the denominator is simply due to 
the I times larger interval of integration in the definition 
of 7c . If H(0, 0) has inversion symmetry, iJ(<E>f, <I>|) is 
transformed to #(—$4, —$4) under inversion. As a con- 
sequence 7c = — 7c (modulo 2tt). This means that in 
systems with inversion symmetry 7c = or ir (modulo 
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2ir) as in previous cases [p]]3],|^,|l(| . 

The values of 7 C and 7 S are easy to predict for systems 
and parameters in which the particles are localized or 
the relevant kinetic terms are not affected by the flux. 
Examples are the CDW and SDW phases with maximum 
order parameter at half filling, for which 7 C = 7s = 
and 7 C = 7 S — n respectively |l0|. Another examples 
which are of interest because of its competition with S 
states, are states with phase segregation (PS) in which 
the particles group together. Let us take a chain of L = 
AN a sites in which the first N a ones are doubly occupied 
and the other are empty (n/l = 1/2). When Eq. (1) is 
applied to this state, all factors give 1, except the last 
one which determines the phase: 7 C = ^2f=i8irj/L = 
7r(7V CT + 1) = (tt) if N a is odd (even). The result is the 
same for states translated to any other place in the chain, 
and to linear combinations of these states (in particular 
with well defined K). This result suggests that choosing 
appropriately L one might be able to detect the boundary 
between S and PS states. Another example in which 
7 C is easy to calculate is a CDW for even L in which 
every second state is singly occupied. Eq. (1) gives 7 C = 
T,7=i*™/ L = n{L + 2) = 0. The same result is obtained 
if every fourth site is doubly occupied. 

The fact that a Berry phase can jump sharply at the 
boundary between an S phase and a CDW or PS state, 
can be useful in finite-size diagonalizations, where cor- 
relation functions vary smoothly at the transition, and 
traditional calculations, like those based on bosoniza- 
tion and conformal- field theory results [ fL5| , which use 
the compressibility obtained numerically from the en- 
ergy necessary to add and remove two particles, may 
have large finite-size effects. For example, in the phase 
diagram of the ID generalized t — J model including a 
three-site term, a superconducting bubble is predicted 
inside the PS region, which is an artifact of finite-size 
effects in the compressibility fltl . Instead, 7 C does not 
depend on the properties of the system for other densities 
and has a unique value on the PS state. In other words, 
although a real phase transition for most systems of inter- 
est can only take place in the thermodynamic limit, un- 
der certain circumstances, both thermodynamic phases 
are characterized by different topological quantum num- 
bers, accessible in finite-size systems, and with small size 
dependence. This was shown in Ref. jl0| for a half filled 
system. Here I calculate 7 C for the following ID model 
at quarter filling (n/l = 1/2): 



7c(PS)= 0) is the convenient choice. Thus, I took L = 12, 
Nf = Ni = 3, for the numerical calculations. As ex- 
pected from the above discussion, I obtain in the region 
V > \U\ 3> t, where charge-charge correlations dominate 
at large distances 7 C (CDW)= 0. 
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I restrict most of the study to U < 0. To distinguish 
the PS state from the one with dominating singlet super- 
conducting correlations at large distances (S), it turns 
out that the odd N a case mentioned above (leading to 




FIG. 2. Charge Berry phase as a function of V for U = —4 
(dashed line) and U = —1 (dotted line). 

In the S region, 7 C = or 7r depending on the param- 
eters. Fig. 2 illustrates 7 C as a function of V. There 
are several jumps in 7 C which do not provide informa- 
tion on phase boundaries. In particular, there is always 
a jump at V = 0, and for large \U\ there is another one 
near V = V} = -{t 2 /\U\ - 4t 4 /|?7| 3 ). To illustrate this 
jump, consider the case \U\ 3> t,V. In this limit, it 
is known that H reduces to a model of N„ hard-core 
bosons, representing doubly occupied sites, which can 
hop to NN positions with matrix elements t'e ± ^*T+*i) 5 
with t' = 2t 2 /\U\ (t' — —2Vc if fourth order corrections 
are included), and NN repulsion V = 2(2V + t'). In 
ID, each hard-core boson can be mapped into a spinless 
fermion with an appropriate change in the boundary con- 
ditions [^] and this model is equivalent to an XXZ model 
Z)<y> a JaSfSf with N a spins up, and J x = J y = 2t' 
, J z = V' . Another way to reach the same result is 
to perform first the unitary transformation = Cjf, 
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(— l)^ct| (which in H has the effect of changing 



the sign of U and $| and transforming the NN repulsion 
in 4V V" , - ^ SfS?;), and then eliminate terms linear in 
t through a canonical transformation. For V' = the 
model can be solved trivially since the interaction van- 
ishes in the spinless-fermion model. Then, one might 
draw a picture like Fig. 1 to represent |<7o($, ^)) an d 
again there is a crossing of levels at $ = ±27r of states 
with the same K = 0, but different (singly) occupied (ef- 
fective) one-particle states k, leading to an undefined j c . 
However, in contrast to Fig. 1, states with K = ±7r/2 
appear in the low-energy manifold and the ground-state 
energy as a function of flux $ displays periodicity in 7r 
(with minima at $ = 7r x integer): the anomalous flux 
quantization characteristic of the S phase |l(],[l^,[l1| . A 
small V leads to a well defined 7 C which depends on the 
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sign of V . I also find other jumps of 7 C for V = and 
at other points inside the S phase, which do not provide 
information on phase boundaries. The physical meaning 
of the other jumps inside the S phase is not obvious. 
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FIG. 3. Phase diagram of the extended Hubbard model 
Eq. (2) determined through jumps in y c , showing the re- 
gions of dominating superconducting (S) or charge-density 
wave (CDW) correlations at large distance, and the region of 
phase segregation (PS). 

Instead, I find that for \U\ 3> t, the first and the 
last transition as a function of V agree very well with 
the limits of the S phase, obtained in the thermody- 
namic limit through the exact Bethe-ansatz solution of 
equivalent XXZ models: from the solution of the above 
mentioned XXZ model fl9[] , we know that PS occurs for 
J z < —J x (ferromagnetic Ising-like case), what in terms 
of the parameters of Eq. (||) means V < Vps — —2t 2 /\U\. 
On the other hand, when V ^> t 2 /\U\, extending previ- 
ous ideas [|l7],[l8| , the entity composed of a doubly occu- 
pied site and an empty site right to it can be mapped 
into a spinless fermion. After an appropriate canoni- 
cal transformation jil| , I find that the resulting spinless 
fermion model has on-site energy U — At 2 /(V — U), NN 
repulsion V = U 2 [l/{V - U) - 1/(3 V - U)}, hopping 
i'e^^T+fi), w ith t' = 2t 2 /(V - U), number of particles 
N' = N a = N/2 and number of sites V = L - N/2. I 
have calculated the correlation exponent K p of this model 
in the thermodynamic limit from its compressibility and 
charge velocity v c mapped appropriately to the original 
model @. The energy and v c were obtained solving the 
corresponding integral equations of the equivalent XXZ 
model jL9| . I obtain that the boundary K p = 1 between 
dominating S or CDW correlations at large distances cor- 
responds to J z = 0.25J X , leading to Vcdw — 0.20\U\. 
For 5< < \U\ < lOt, I find that for L = 12, the first jump 
in 7 C with increasing V is above Vps by ~ 25% and the 
last jump is above Vcdw by ~ 15%. These numbers are 
reduced by a factor ~ 2/3 if fourth order corrections to t' 
and V are included. I believe that the remaining discrep- 
ancy is due to finite size effects. For smaller values of \U\, 
the first and last jumps in 7 C , shown in Fig. 3 agree with 



previous numerical calculations of K p and compressibil- 
ity |l8| , |2"o| ] . In correspondence with these calculations, I 

find an island for high V and small \U\ (V ~ 8t, U 2t, 

not shown) in which j c jumps to 7r and K p > 1. How- 
ever, using again an appropriate mapping to two XXZ 
models, I obtain that this is a finite-size effect and that in 
the thermodynamic limit the system phase separates into 
two phases with dominating CDW correlations, in agree- 
ment with previous suggestions |l8). One phase is the 
one just described for V > Vcdw, and in the other one 
fermion objects composed of one particle and an empty 
site next to it move with hopping t and NN attraction 
V = -2t 2 /V. 

Another limit in which the charge dynamics is equiv- 
alent to that of interacting spinless fcrmions is that of 
U = +oo. Following similar methods as those already 
used for V = |pl| , the model can be mapped into the 
spin polarized case iVj = L/2, N^ = 0, with a shift 
A$ in the BC. Calculating K p in the thermodynamic 
limit, I find that superconducting correlations dominate 
for — 2t < V < Vcdw, where up to three digits accuracy 
Vcdw = -a/2*- For 8 < L < 16, 7c jumps exactly at 
this value . Instead, 7 C does not detect the transition to 
the PS state at V = —2t in this case. 

In the rest of this Letter, I discuss the spin Berry phase 
7 S . As noted earlier |l^,|ll| , 7 C is transformed into 7 S + 7r 
(of a different Hamiltonian in general) and vice versa by 
the transformation = c^, = ( — l) J ct^ . Thus, the 
definition of 7 S can be generalized to JV-f ^ N± in the 
same way as 7 C , keeping $j = — $| in the I cycles, and 
I should be such that (Nf — N\)l/L is an integer. It is 
also clear from this transformation and what I obtained 
for 7c (Eq.(|| )) that A(P T - P{) = {e/2Trl)A ls (mod- 
ulo e/Z), where P CT is the macroscopic polarization of the 
particles with spin a. If N\ = iVj_, as in most cases of 
interest, the previous definition of "f s with Z = 1 p0| , is 
valid for any filling, and for ID systems with spin SU(2) 
symmetry, 7 S jumps in ir when a spin gap A s opens. To 
show this, I use recent results from continuum limit the- 
ory and bosonization [p2[ , which show that the opening 
of a spin gap is determined by the crossing of the lowest 
energy states of H with total spin S = and S = 1 within 
the K = sector, for periodic (antiperiodic) BC if N a is 
even (odd). These BC correspond to a particular point 
<t> = <t> c with <t> c = ($ c = tt) in the trajectory with 
K = K = of H ($, — $) used to define j s . Usually this 
point corresponds to the maximum of the ground-state 
energy Eq(Q, — $) as a function of while its minimum 
lies at $ = 7r (<f> = 0) for N a even (odd), as expected 
from the non- interacting limit. 

From its definition, it is obvious that to calculate 7 S 
any (single valued) ^-dependent gauge can be chosen. I 
take the representation in terms of the operators c -^ = 
e ij* c /L c t ^ gt _ g-jj^c/igt m sucn a wa y tnat m 

resulting Hamiltonian H($>, — $), all hoppings are real for 
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<f> = $ c (as in H"(<I> C , — $ c ) 
can approximate — $) 



Sufficiently near to <E> C , one 
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where the derivative is evaluated at $ c . While at $ = $ c , 
the SU(2) symmetry is conserved, for <I> ^ <E> C only spin 
rotations around the quantization axis z are conserved 
by H. Specifically, dH/d<& is equal to the paramagnetic 
part of the difference between up and down currents, and 
changes sign under spin rotations in 7r around the x or y 
axis. I denote these rotations by R and call \s), E s ( \t), 
Et) the lowest singlet (triplet with spin projection zero) 
eigenstate of _ff($ c , — $ c ) and its energy. It is known 
that R\s) = \s), and R\t) = -\t). Since dH/d® is odd 
under R, this implies (s\dH /d<f>\s) = (t\dH/d$\t) = 0, 
However, in general (s\dH / d<f>\t) = A ^ 0, where A can 
be made real by a suitable change of the phase of one of 
the states. Near the crossing point between lowest singlet 
and triplet energies, the ground state |<?o($> — $)) results 
from the diagonalization of the low-energy part Hle of 
H containing the states \s) and \t), which according to 
the above symmetry arguments, has the form: 



H LE (Q,-$) 



f E s 



(4) 



It is easy to see that this leads to a jump in n in "f s 
when Et crosses Es- For example, one can exploit again 
the gauge invariance to render (s\go(§>, — $)) real and 
positive for $ near $ c , leading to \go($ c + + ,— $ c — 
0+)) = | 5o ($ c - 0+,-$ c + 0+)> for E s < E T but 
| ff0 (* c + 0+,-$ c - 0+)) = -|flo(*c - 0+,-$ c + 0+)) 
for Es > Et and no other contribution to 7 S in the 
vicinity of <i> c . I obtain numerically in the extended Hub- 
bard model with or without correlated hopping ||io|| that 
7 S = in the spin gapped phase and therefore j s = tt if 
A s = 0. According to continuum limit theory and renor- 
malization group, in ID, when K p > 1, the opening of 
a spin gap signals the transition from dominating triplet 
to singlet superconducting correlations at large distance. 

In summary, I have shown the ability of two Berry 
phases, and topological quantum numbers derived from 
them, to detect phase transitions, particularly in strongly 
correlated systems with non-integer number of particles 
per site, for which no previous applications of Berry 
phases were considered. I hope that this study will stimu- 
late further research and applications of the Berry phases 
in phase transitions, particularly in dimensions higher 
than one, where alternative methods which use results of 
conformal field theory, or exact solutions p^ , p^ , p^| - p0|j22| 
are not available. 

I am indebted to Gerardo Ortiz for providing me with 
details of his calculations in Ref. || and useful discus- 
sions. The numerical diagonalization was done using 
subroutines developed by Eduardo Gagliano, who un- 
fortunately died on February 12, 1998. I am partially 
supported by CONICET, Argentina. 
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